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CN . ABSTRACT 

We present a cross-correlation of the European Large Area ISO survey (ELAIS) with 
^ ! the ELAIS Deep X-ray Survey of the Nl and N2 fields. There are 7 Chandra point 

sources with matches in the ELAIS Final Analysis 15/ini catalogue, out of a total of 
28 extragalactic ISO sources present in the Chandra fields. Five of these are consistent 
with AGN giving an AGN fraction of ~ 19 per cent in the 15 nm flux range 0.8 — 6 
i-S^ ' mjy. We have co-added the hard X-ray fluxes of the individually-undetected ISO 

sources and find a low significance detection consistent with star formation in the 
Q . remaining population. We combine our point source cross-correlation fraction with 

$H ' the XMM-Newton observations of the Lockman Hole and Chandra observations of the 

, Hubble Deep Field North to constrain source count models of the mid-infrared galaxy 

■ population. The low dust-enshrouded AGN fraction in ELAIS implied by the number 

of cross-identifications between the ELAIS mid-infrared sample and the Chandra point 
sources is encouraging for the use of mid-infrared surveys to constrain the cosmic 
' star formation history, provided there arc not further large undetected populations of 

H ' Compton-thick AGN. 

! 

Key words: surveys - X-rays: general - X-rays: galaxies - galaxies: active - quasars: 
general 
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1 INTRODUCTION 

Enormous progress has recently been made in resolving the 
sources that comprise the extragalactic hard X-ray back- 
ground (e.g. Mushotzky et al. 2000, Cowie et al. 2002, 
Moretti et al. 2003). Studies are now focusing on charac- 
terising the nature of these sources through multi- waveband 
imaging and spectroscopy. Surveys with the ISO satellite 
have recently discovered strong evolution in the mid-infrared 
galaxy population (e.g. Elbaz et al. 1999, Serjeant et al. 
2000, Chary & Elbaz 2001, Gruppioni et al. 2002) which is 
thought to be mainly due to a strongly evolving obscured 
cosmic star formation history (e.g. Aussel et al. 1999), al- 



though there is also a contribution from dust-enshrouded 
AGN. Optical spectroscopic follow-ups of these samples are 
underway (e.g. La Franca et al. 2003, Gonzalez- Solares et al. 
2004, Perez-Fournon et al. in preparation) but as the sources 
might not be optically-thin at optical wavelengths, nor emit- 
ting isotropically, it is possible that heavily dust-enshrouded 
populations may be mis-classified by this approach. The 
Chandra observations of the HDF-North (e.g. Hornschmeier 
et al. 2001, Alexander et al. 2002) placed constraints on 
the fraction of dust-enshrouded AGN at the faintest end of 
the 15/im source counts. Most of the known 15^m galax- 
ies are at somewhat higher flux densities, however. Some 
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inroads were made on the AGN fraction for brighter mid- 
infrared sources by Alexander et al. (2001) using BeppoSAX 
observations of the ELAIS survey (described below), but the 
BeppoSAX depth was not sufficient to detect most Compton- 
thin Seyfert II galaxies in the targeted ELAIS field. Fadda 
et al. (2002) provided the first study with a reasonably sig- 
nificant number of sources with hard X-ray and mid-infrared 
emission. They combined XMM-Newton observations of the 
Lockman Hole and Chandra observations of the HDF-N with 
coincident ISOCAM data. This paper will provide further 
statistics using 2 Chandra pointings in the ELAIS north- 
ern survey regions Nl and N2. An upcoming analysis of the 
Spitzer Wide-area Infrared Extragalactic Survey (SWIRE, 
Lonsdale et al. 2003) observations in the Nl region will pro- 
vide enhanced statistics in complementary wavebands. 

The European Large Area ISO Survey, ELAIS, was the 
largest open time project on ISO, covering wavelengths from 
6.7/im to 175/im (Oliver et al. 2000, Rowan-Robinson et al. 
2004). The survey, and its follow-ups, have many ambitious 
aims, including tracing the cosmic star formation history to 
z ^ 1 and the discovery of ultra- and hyper-luminous galax- 
ies at high redshift. The mid-infrared source counts (Serjeant 
et al. 2000, Gruppioni et al. 2002, Elbaz et al. in preparation) 
and far-infrared source counts (Efstathiou et al. 2000) both 
show evidence for strong evolution, as does the far-infrared 
luminosity function (Serjeant et al. 2001). The ELAIS mid- 
infrared source counts cover the transition from Euclidean 
slope to steep evolution; by virtue of the large area and large 
investment in observing time, ELAIS represents by far the 
largest sample of galaxies from this strongly evolving mid- 
infrared population. A high proportion have been revealed 
as ultra- luminous infrared galaxies (14% of 15/im galaxies 
with known z) including 9 hyper-luminous infrared galaxies 
(Morel et al. 2001, Rowan-Robinson et al. 2004). 

Thanks to extensive multi-wavelength coverage, the 
ELAIS fields have now arguably become the most well stud- 
ied regions of their size, and natural targets for on-going or 
planned large-area surveys with the most powerful ground 
and space-based facilities. Further details on ELAIS multi- 
wavelength observations and catalogues are presented in 
Rowan- Robinson et al. (2004). In particular, ELAIS 15 
observations will complement the SWIRE survey in three ar- 
eas (Nl, N2 and SI) by covering the 8-24um gap in Spitzer^s 
filters. 

In the ELAIS Deep X-ray survey we made deep Chan- 
dra pointings of the ELAIS Nl and N2 fields to a limiting 
depth of ~ 10"^^ erg s"^ cm"^ in the 0.5 - 8keV band, 
with a total exposure of ~ 75ks in each field. Manners et al. 
(2003) present the data and source counts. Gonzalez-Solares 
et al. (in preparation) present follow-up imaging and optical 
spectroscopy of the sample and Willott et al. (2003) present 
Subaru infrared spectroscopy. Almaini et al. (2003) present 
a cross-correlation of the Chandra sources with the sub-mm 
sources of Scott et al. (2002) and Fox et al. (2002), and 
measure the clustering of the Chandra population. 

Here we present a cross-correlation of the Chandra X- 
ray sources with the mid-infrared sources from the ELAIS 
survey. Section |21 summarises the data acquisition in the 
mid-infrared and X-ray, and provides references to more ex- 
haustive descriptions for the interested reader. Section|H]de- 
scribes the cross-correlation between the ISO and Chandra 
data, and section 2] discusses the significance of the results. 
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Figure 1. Positions of the 28 extragalactic 15 /im ISO sources 
(black circles) superimposed on the smoothed 0.5 - 8 keV Chan- 
dra images. The images have been adaptively smoothed using the 
flux-conserving algorithm CSMOOTH from the Chandra Interac- 
tive Analysis of Observations (CIAO) package. 



2 OBSERVATIONS 

As part of ELAIS, ISOCAM observations were carried out 
covermg 10.3 deg^ at 15 ^m. The observations were carried 
out in raster mode, with most of the survey area confined 
to three northern fields (Nl: 2.67 deg^, N2: 2.67 deg^ N3: 
0.88 deg^) and one southern field (SI: 3.96 deg^). Further 
details on the observation strategy can be found in Oliver et 
al. (2000). 

Data reduction of ELAIS 15 ^m observations was re- 
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Table 1. Chandra sources with ISOCAM 15 fim counterparts in the ELAIS regions Nl & N2. 1 Chandra source ID from Manners et 
al. (2003). ^Manners et al. (2003). ^Rowan-Robinson et al. (2004). ^Offset between Chandra and 15 /im ISO positions. ^Probability of 
this cross-correlation being a random association (see section 0. ^Chandra full band flux (Manners et al. 2003). ^ISOCAM 15 /xm flux 
density (Rowan-Robinson et al. 2004). * Spectroscopic redshifts (t photometric, Icr error ~ 0.1) reported in Rowan-Robinson et al. (2003). 
^Likely object class as suggested by this paper ("section 14. li or -fpreviously confirmed by optical spectroscopy. 
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Table 2. Multi-waveband data for each cross-correlated source. Uncertainties are Icr. X-ray data are from Manners et al. (2003). Other 
waveband data are taken from Rowan-Robinson et al. (2004) except ^^Moshir et al. (1990), 2Gonzalez-Solares et al. (in preparation), 
^Willott et al. (2003), ''Ciliegi et al. (1999), ^aperture magnitudes provided by E. Gonzalez-Solares. Also included are the mid-IR to 
hard X-ray relation («/x) ^'^<^ the X-ray hardness ratio (HR) as defined in Manners et al. (2003). Units are "mJy, ''Vega magnitudes, 
and ''lO^i^ erg cm~2 g-i 



cently completed by Vaccari et al. (in preparation). Data 
reduction was carried out using the LARI method (Lari 
et al. 2001, Lari et al. 2003), a new technique devised for 
the reduction of ISOCAM and ISOPHOT imaging data. 
Based on a physical model of the ISO detectors' behaviour, 
the method is particularly suited for the reliable detection 
of faint sources in ISO surveys, allowing sensitivity to be 
pushed to the instrumental limits. 

A sample of 1056 sources (490 in Nl and 566 in N2) 
were detected with S/N > 5, spanning the 0.5 - 100 mJy 
flux range and thus filling the gap between the ISOCAM 
deep surveys (e.g. Elbaz et al. 1999) and the IRAS Faint 
Source Catalogue (Moshir et al. 1990). 

The ELAIS Deep X-ray Survey was centred on regions 
in the ELAIS Nl and N2 fields selected to have low cirrus 
and HI column density (see Oliver et al. 2000 for more de- 
tails). The Chandra observations are described in detail else- 
where (Manners et al. 2003), though we summarise the main 



points here. The observations were taken with the Chandra 
ACIS (Advanced CCD Imaging Spectrometer) array. Inte- 
grations of 75ks were taken in each of Nl and N2. The 
pointing centroids are 16:10:20.11 +54:33:22.3 in Nl and 
16:36:46.99 +41:01:33.7 in N2, with an area of 16.9' x 16.9' 
covered by the ACIS-I chips in each case, giving a total area 
of 571 sq. arcmin. The limiting flux levels are 4.6 x 10"^'' 
erg s"i cm~^ in the 0.5 - 2 keV band, and 2.2 x 10"^^ erg 
s"'^ cm~'^ in the 2-8 keV band. The Nl region contains 125 
Chandra sources from the 4 ACIS-I chips and 5 from the 
ACIS-S2 chip. The N2 region has 99 sources on the ACIS- 
I chips and 4 on the ACIS-S2 chip. Only sources from the 
ACIS-I chips were used for this analysis due to the poor 
resolution of the off-axis ACIS-S chips. 
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Figure 2. Chandra sources with a 15 /im counterpart in tlie Nl region. 25" X 25" r'-band postage stamps are from a 5400 second 
exposure witli the INT WFC and are displayed in positive and negative grey-scale for clarity. Black or white circles indicate position 
of the ISO source and grey ellipses give the 3fT size of the X-ray PSF. Spectral energy distributions are at rest-frame wavelengths and 
displayed in units of luminosity where redshifts are available. 



3 POSITIONAL CORRELATION BETWEEN 

THE CHANDRA AND ISO SOURCES 

In the Nl region there are 9 15/im Final Analysis ELAIS 
sources within the Chandra region, three of which are iden- 
tified with stars and the remainder of which are galaxy IDs. 
In the N2 region there are 24 15/im sources, two of which 
are identified with stars and the rest with galaxies or blank 
fields in the r' band image. The positions of the 28 extra- 
galactic sources are plotted in Fig. superimposed on the 
smoothed Chandra images. 

We performed a simple near-neighbour search to cross- 
correlate the extragalactic Chandra and ISO sources within 
the area of the Chandra ACIS-I chips, using a 5" search 
radius. Astrometric la errors for the Chandra sources are 
~ 1" (Manners et al. 2003), while the nominal astrometric 



accuracy for the ELAIS 15 /tm FA catalogue sources ranges 
from ~ 0.8" to 2.0" (Vaccari et al. in preparation). The 
search radius of 5" was chosen as the approximate sum of 
the 2(T astrometric errors. Three matches were found in the 
Nl region and four in the N2 region (table Q Fig. 0. All 
seven matches are with high reliability ISO sources (5cr). 

To ensure the associations were real we calculated the 
probability of a random association between each mid-IR 
source and its X-ray counterpart. Following Fadda et al. 
(2002), we assume the X-ray counterpart belongs to a Pois- 
sonian distributed population so that 

P = l_e-N(>S)-rf' (1) 

where P is the probability of a random association within 
an offset distance d. N(>S) is the number density of sources 
with flux greater than the possible X-ray counterpart (S). 
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Figure 3. Chandra sources with a 15 //m counterpart in the N2 region. 25" x 25" r'-band postage stamps are from a 7800 second 
exposure with the WHT and are displayed in positive and negative grey-scale for clarity. Black or white circles indicate position of the 
ISO source and grey ellipses give the 3cr size of the X-ray PSF. Spectral energy distributions are at rest-frame wavelengths and displayed 
in units of luminosity. 
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We calculate this probability for each source (table with 
reference to the log(N)-log(S) relation for these regions 
(Manners et al. 2003). The chances of random associations 
are found to be very low. 

Figs. |5| & 13 display the cross-correlations overlaid on 
optical postage stamps. 



4 RESULTS AND DISCUSSION 
4.1 Identifications 

Properties of the 7 cross-correlated sources are reported in 
table Multi-waveband data are reported and referenced 
in table 121 Also reported in table |21 are the mid-infrared to 
X-ray spectral indices [aix, described in section [4.21 and 
the X-ray hardness ratios (HR, from Manners et al. 2003). 
r' band postage stamps to a depth of r' ~26 are shown in 
Figs. |5| & 13 together with the spectral energy distribution 
(SED) of each source. For the 6 sources with available red- 
shifts, SEDs are displayed in the rest frame in units of lumi- 
nosity assuming a cosmology with Q,a = 0.73, Q,m — 0.27, 
and Ho = 71 km s~^ Mpc~^. For the source without a red- 
shift (Nl_49), the SED is displayed in the observed frame in 
units of flux. In order to derive flux densities for the SEDs, 
magnitudes were converted using the following zero points 
for each photometric band: K 657 Jy; H 1020 Jy; J 1600 Jy; 
i' 2491 Jy; r' 3133 Jy; g' 3876 Jy; U 1810 Jy. Flux densities 
at 1 keV and 5 keV were derived from the 0.5-2 keV and 
2-8 keV Chandra bands respectively, assuming a power- 
law spectrum with photon index F = 1.7 within each band. 
A description follows of the properties and most probable 
identifications for each source: 

Nl_20 (CXOENl J161046. 5+543538) A complex r' 
band morphology indicative of a recent merger. This low 
redshift {z = 0.0634) object has a soft X-ray spectrum and 
is below the detection threshold in the 2-8 keV band im- 
age. The mid-infrared to X-ray spectral index is also quite 
steep {aix ~ 1.51). These characteristics are consistent with 
identification as a starburst galaxy (SB, table 0. 

Nl_49 (CXOENl J161020. 8+543900) Spectroscopic 
identification of this object is hampered by the presence of 
a nearby bright star (as can be seen in the r' band postage 
stamp). No redshift is available, however a very hard X-ray 
spectrum combined with a relatively flat mid-infrared to 
X-ray spectral index indicates likely identification with an 
AGN. The X-ray spectrum is consistent with an apparent 
absorbing column of Nh ~ 6 ± 2 x 10^^ cm~^ (assuming 
an underlying power law of a = 0.7 and z = 0). Depending 
on the redshift of this source the actual Nh is likely to be 
higher with apparent absorbing column scaling as {1 + z)'^'^. 
The mid-infrared to X-ray spectral index {aix ~ 1.32) is 
also consistent with an absorbed AGN (see section 14.211 . 
Given an AGN 2 classification in table Q 

Nl_69 (CXOENl J161003. 1+543628) Spectroscopic 
redshift of 0.2675. The flat mid-infrared to X-ray spectral 
index, luminosity and X-ray hardness ratio are all typical 
of a Seyfert galaxy (AGN 1, table [TJ. 

N2_25 (CXOEN2 J163655. 7+405910) 
Spectroscopically confirmed quasar at a redshift of 2.61 
(AGN 1, table 0. Its properties are extensively reported in 
Willott et al. (2003). 
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Figure 4. Mid-infrared to hard X-ray spectral indices (aix) 
vs. redshift. Shaded areas represent the error boxes for template 
SEDs of QSOs, Seyfert 2s and Starburst galaxies as reported in 
Alexander et al. (2001, Fig. 7). cf. Fadda et al. (2002, Fig. 8). 

N2_33 (CXOEN2 J163651. 6+405600) 

Spectroscopically identified as a Seyfert 2 galaxy at a 
redshift of 0.4762 (AGN 2, table[TJ. The r' band image indi- 
cates an interaction is occurring with a smaller neighbour. 
Radio morphology suggests this to be an FR-II type AGN. 

N2_52 (CXOEN2 J163629. 7+410222) Elliptical 
galaxy with no obvious point-like nuclear source. A pho- 
tometric redshift of ~ 0.02 suggests an X-ray luminosity 
that is low for an active galaxy (1 x 10*" erg s~^ at 1 keV, 
2 X 10*" erg s~^ at 5 keV). However, mid-infrared to X-ray 
spectral index and X-ray hardness ratio are both consistent 
with an AGN (AGN 1, table Q. 

N2_107 (CXOEN2 J163608. 4+410507) Interesting r' 
band morphology displaying merging galaxies with a double 
nucleus. Spectroscopic redshift of 0.1683. The X-ray spec- 
trum is soft and of low-luminosity (6 x 10*" erg s~^ at 1 
keV) with the source undetected in the 2-8 keV band im- 
age. The relatively luminous infra-red spectrum supports 
identification with a starburst galaxy (SB, table 0. 

Of the 7 matched sources, 1 is a spectroscopically con- 
firmed quasar, 1 is spectroscopically identified as a Seyfert 2, 
a further 3 display properties of AGN, and 2 have properties 
consistent with starburst galaxies. 

4.2 The mid-infrared to X-ray spectral indices 

Values for the mid-infrared to X-ray spectral index (aix) 
are calculated using the fiux density observed at 15/im and 
5 keV, assuming a power-law spectrum of the form F^ oc 
p-"ix fpjjg density at 5 keV is calculated from fiux in 
the 2-8 keV Chandra band. For the 2 objects undetected 
in this band, the full 0.5 - 8 keV band is used. 

The mid-infrared to X-ray spectral index can be a useful 
indicator to distinguish between AGN and starburst galax- 
ies (e.g. Alexander et al. 2001, Fadda et al. 2002). Starburst 
galaxies are found to have high values of aix, while type-1 
AGN have low values {ctix < 1.2). Type-2 AGN have val- 
ues in between depending on the amount of obscuration. 
Fig. El plots the aix values for our 6 matched sources with 
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available redshifts. The figure displays values of aix as a 
function of redshift for template SEDs compiled by Alexan- 
der et al. (2001). Source N2_33, spectroscopically identified 
as a Seyfert 2 lies squarely within the region of the Seyfert 2 
templates. N2_25, spectroscopically identified as a QSO has 
a value slightly higher than the QSO templates although 
still lower than expected for a Seyfert 2. Of the remain- 
ing objects studied here, Nl_69 and N2_52 lie in the region 
of type-1 AGN, Nl_20 and N2_107 are consistent with star- 
bursts or highly obscured AGN, and Nl_49 is consistent with 
a type-2 AGN. 

4.3 Constraints on models of mid-infrared source 
counts 

In the Nl region there are 6 extragalactic 15/im Final Anal- 
ysis ELAIS sources within the Chandra region. In the N2 
region conversely there are 22 extragalactic 15/im sources. 
The Chandra N2 region falls within an area of repeated ISO 
observations while the Chandra Nl region does not. The 
difference in source counts, however, is larger than expected 
and may indicate clustering on scales larger than the field 
size. This highlights the need for larger areas to be surveyed 
before sufficient count statistics can be gained. Here we con- 
sider our results along with number counts from the Hubble 
Deep Field - North (HDF-N) and the Lockman Hole (Fadda 
et al. 2002) in order to provide a comparison with recent 
models of mid-infrared source counts. 

The 15pim flux range 0.8 — 6.0 mjy is well covered by the 
ISO observations in our Chandra regions. In this fiux range 
there are 26 extragalactic 15/im sources in our sample. We 
take 5 of these objects to contain evidence of AGN from 
the Chandra and multi-waveband data. This gives an AGN 
fraction of 5/26 extragalactic sources, or 0.19 ± 0.09 over 
the given flux range, in a total area of ~ 571 sq. arcmin. 
This fraction will be a lower limit if any heavily obscured, 
Compton-thick AGN are present in our sample. 

Fadda et al. (2002) report AGN number counts for mid- 
infrared sources in the Lockman Hole and HDF-N. They find 
for the Lockman Hole an AGN fraction of 13/103 extragalac- 
tic sources (0.13±0.04) within the 15/im fiux range 0.5 — 3.0 
mJy in a total area of 218 sq. arcmin. For the HDF-N the 
AGN fraction is 5/42 (0.12 ± 0.05) over a 15/im fiux range 
of 0.1 — 0.5 mJy in a total area of 24.3 sq. arcmin. 

Wc compare these results with models for the AGN con- 
tribution to mid-infrared surveys by Pearson (2001) and 
King & Rowan- Robinson (2003). These models primarily 
attempt to model the star formation history by fitting in- 
frared source count observations with contributions from 
normal galaxies, starbursts, ultra-luminous infrared galax- 
ies (ULIRGs) and AGN. The Pearson (2001) model uses 
local luminosity functions and pure luminosity evolution to 
describe the normal galaxy, starburst and AGN populations 
while the ULIRG component is evolved in both density and 
luminosity. The 15/im and 850/im source counts are used to 
constrain the model. The best fit to the observations is ob- 
tained where the ULIRG component undergoes two major 
phases of evolution, rapid merging to z ~ 1 and an expo- 
nential evolution in luminosity to higher redshifts. This has 
an effect on the AGN fraction of the source counts which 
will be tested here. The King & Rowan-Robinson (2003) 
models are updated versions of models developed by Rowan- 
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Figure 5. AGN fraction in the mid-infrared population derived 
from the ELAIS Deep X-ray survey (this work), as well as frac- 
tions derived by Fadda et al. (2002) from Chandra observations 
of the Hubble Deep Field North (HDF-N) and XMM-Newton ob- 
servations of the Lockman Hole (LH). Also plotted are the predic- 
tions from the integral counts of Pearson (2001), and two models 
from King & Rowan-Robinson (2003) (based on the models of 
Rowan-Robinson 2001). 



A 



o 




Pearson 2001 
Manners 2002 
King 2003 (Lambda model) 
King 2003 (Einstein de Sitter model) 



3L 
0.01 



0.10 



1.00 



10.00 



(mJy) 



Figure 6. Integral AGN number counts for 15/tm surveys. Values 
for the Lockman Hole and HDF-N are taken from Fadda et al. 
(2002). Number count predictions are taken from the same models 
displayed in Fig.|^ plus that of Manners (2002). 



Robinson (2001) and are calculated for lambda and Einstein 
de Sitter cosmologies. They allow for both density and lu- 
minosity evolution in all four populations. The AGN con- 
tribution in both models is found using a 12/im luminosity 
function from Rush, Malkan & Spinoglio (1993) which does 
not include Compton-thick AGN and can therefore be di- 
rectly compared with these observations. Fig.l^displays the 
AGN fraction predicted by these models against source fiux. 
All 3 models shown are consistent with AGN fractions re- 
ported here for the ELAIS regions and for the HDF-N. The 
AGN fraction in the Lockman Hole is slightly lower than the 
predictions of these models. 

We also compare the observations to modelled pre- 
dictions for the cumulative AGN source counts at 15/im 
(Fig. |SJ. To make a reliable comparison it was necessary 
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to correct our AGN source counts for completeness. While 
the 15 fira data reduction has been completed by Vaccari et 
al. (in preparation), completeness estimates are to be pub- 
lished in Lari et al. (in preparation) and are not yet available. 
Therefore, values from Gruppioni et al. (2002) were taken as 
representative. By estimating the completeness at the flux 
level of our sources we calculate the expected number of 
AGN above the flux limit to be ~ 7.7 (ELAIS data point, 
Fig.|SJl. We can now compare a further model from Manners 
(2002) which assumes AGN are responsible for the entire 
hard X-ray background and predicts the emission of these 
AGN in the mid-infrared. Observations in the HDF-N agree 
with models by King & Rowan- Robinson (2003) and Man- 
ners (2002), while the Pearson (2001) model over-predicts 
the number of sources at this depth. AGN counts from the 
Lockman Hole are consistent with the Pearson (2001) and 
King & Rowan-Robinson (2003, Lambda) models, whilst the 
AGN counts reported here for the ELAIS survey agree well 
with all the model predictions. 

4.4 Statistical hard X-ray limits of ELAIS ISO 
sources 

One further method of constraining the source count mod- 
els is the statistical detection of ISO sources in the Chan- 
dra map. By co-adding the Chandra fluxes at the positions 
of ISO sources not detected individually, we can obtain a 
constraint on the mean hard X-ray flux and hard X-ray to 
mid-infrared flux ratio for the remaining population. 

In the Nl region there are 3 extragalactic sources un- 
detected by Chandra whilst in the N2 region there are 18 
(see Fig. 0. In order to obtain limits on the X-ray flux for 
these sources it was necessary to define source regions on the 
Chandra images. These regions were positioned accurately 
using the astrometry obtained by cross-correlating Chandra 
and r'-band sources (Manners et al. 2003). The size of each 
region was defined as a circle containing 95% encircled en- 
ergy for a monoenergetic Chandra PSF at 4.51 keV at the 
relevant source position. 

For each field, around 6 circular background regions of 
radii ranging from ~ 1' — 2' were selected. These were cho- 
sen carefully to avoid contamination from known Chandra 
sources and to cover regions of typical effective exposure. 
The background counts expected in each source region were 
then calculated by correcting for the difference in effective 
exposure. This method was preferred over the use of indi- 
vidual local background estimates due to the slow change in 
background counts over the image and the better statistics 
gained from using the larger background regions. 

Counts were extracted from the hard (2-8 keV) and 
soft (0.5 - 2 keV) band Chandra images. Source regions were 
individually corrected for effective exposure, background 
subtracted, and then co-added to obtain an estimate of the 
total flux. 

In the hard band we obtain a total of 58 counts with 
49.3 background counts expected. The Poisson probability 
of obtaining 58 or more counts with 49.3 expected is 0.123, 
equivalent to a tentative detection at a confldence level of 
87.7%. This is equivalent to a mean flux of 2.2 ± 1.5 xlO"" 
erg cm~^ s"'^ for each of the 21 sources. In the soft band we 
obtain a total of 43 counts with 24.1 background counts 
expected. The Poisson probability of obtaining 43 or more 



counts with 24.1 expected is 3 x 10"'*, equivalent to a reliable 
detection at a confidence level of 99.97%. These counts are 
equivalent to a mean flux of 6.7 ± 2.5 xl0~^^ erg cm"'^ s~^. 

In order to check the validity of the stacking analysis we 
repeated the procedure with the positions randomised over 
the areas of the Chandra field that are free of known sources. 
This was repeated 10 times for each band. The distributions 
of the counts obtained were consistent with a Gaussian of 
mean and variance equal to the expected number of back- 
ground counts in each case. 

The stronger detection in the soft band is indicative of 
the higher efficiency of Chandra in this band and does not 
imply these sources have a particularly soft X-ray spectrum. 
The hardness ratio (defined in Manners et al. 2003) for the 
co-added regions is -0.09 ± 0.38, somewhat harder than ex- 
pected for starbursts. However, the low significance of the 
hard X-ray detection means this must be used with caution. 

Seven of the unmatched sources have spectroscopic red- 
shifts available, ranging from z = 0.10 to 0.24. At a typical 
redshift of 0.2 the mean fluxes given above would equate to 
luminosities of 5.4 x 10^^ at 1 keV and 1.9 x 10"° at 5 keV. 
This is entirely consistent with a starburst origin for the 
X-ray emission. The mean 15/im flux for the 21 unmatched 
sources is 1.489 mjy. This leads to a mean spectral index of 
aix ~ 1.68, also consistent with purely starburst galaxies. 

It should be noted that 'Compton thick' AGN will not 
appear in the X-ray data and constitute an unknown fraction 
of the unmatched sources. 



4.5 The cosmic star formation history from 
mid-infrared flux limited samples 

Of the sources detected in field surveys by the ISO satellite, 
by far the largest fraction and largest number of moderate- 
redshift {z > 0.5) systems are found in mid-infrared ISO- 
CAM surveys, as opposed to the far-infrared (90 — 175^tm) 
surveys conducted by the ISOPHOT instrument. Also, the 
surveys now being performed by the Spitzer satellite, such 
as SWIRE (Lonsdale et al. 2003), will have high-redshift ob- 
jects preferentially detected in the mid-infrared passbands. 
The AGN torus emission peaks in the infrared (e.g. Haas 
et al. 1998) so it might be expected that a large popula- 
tion of AGN would occur in these surveys in addition to 
high-z star forming galaxies. Our AGN fractions, together 
with other studies (Fadda et al. 2002, Alexander et al. 2002), 
show this population to be a fairly well-determined minor- 
ity and are encouraging for the use of mid-infrared samples 
for constraining the cosmic star formation history. The mid- 
infrared luminosity is a reasonably good star formation rate 
indicator, albeit affected by complicated K-correction effects 
(e.g. Xu et al. 1998). 

The main caveat is that our Chandra data still do not 
exclude the possibility of a large population of Compton- 
thick objects at moderate redshifts (2 ~ 1). Alexander et 
al. (2002), however, argue that this fraction should be low. 
They performed an X-ray stacking analysis on those infrared 
galaxies that were not clearly AGN at X-ray energies but 
were still individually detected in the 1 Ms Chandra Deep 
Field North Survey. They find an average X-ray spectral 
slope of r = 2.0, indicating a low contribution from the 
much flatter spectra of obscured AGN. 
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5 CONCLUSIONS 

We performed a cross-correlation of X-ray and mid-infrared 
point sources in the ELAIS Nl and N2 fields, using data from 
the Chandra and ISO satellites. 7 extragalactic matches are 
found (out of a total of 28 ISO sources) within the area 
of the Chandra ACIS-I chips. 2 of these are spectroscopi- 
cally identified as AGN. Based on X-ray to IR flux ratios, 
X-ray hardness ratios, and luminosities, 3 of the remaining 
5 are also consistent with AGN while the other 2 are consis- 
tent with starburst galaxies. This provides an AGN fraction 
of ~ 19 per cent in the 15^m flux range 0.8 — 6 mjy. We 
have co-added the hard X-ray flux at the positions of the 21 
undetected ISO sources providing only a 1.4 sigma detec- 
tion. This translates to a mean hard X-ray to mid-IR flux 
ratio consistent with star formation in these objects. Our 
cross-correlations, when compared with XMM-Newton ob- 
servations of the Lockman Hole and Chandra observations 
of the northern Hubble Deep Field, allow us to place con- 
straints on source count models of the 15^m source popula- 
tion. The AGN fractions and number counts are broadly con- 
sistent with models by Pearson (2001) and King & Rowan- 
Robinson (2003). We argue our data is encouraging for the 
use of mid-infrared samples to constrain the cosmic star for- 
mation history, provided there is not a large contribution 
from Compton-thick AGN. 
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